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GRAVITATIONAL STABILITY OF A MASSLESS SCALAR FIELD 

V. M. Razorenov UDC 530.1(04) 



A solution is obtained to the Einstein equations for a static, spherically sym- 
metric, massless scalar field. The stability of a static, massless scalar field 
for a point source is studied. It is found that the field is unstable for small 
ratios of the scalar charge of the source to its mass but is stable for large 
ratios. It is proven that there exists a dimensionless number which limits the 
region of stability for the field. 



The problem of gravitational stability of a static, spherically symmetric, massless 
scalar field was formulated in [1], where the authors concluded that the scalar field of 
isolated sources is unstable. We will show that this conclusion is unfounded and that there 
exists a region of stability when the ratio of the charge to the mass of a point source ex- 
ceeds some fully specified quantity (there is a dimensionless parameter that enters into 
the static solution which limits the regions of stability and instability). 

We have for the metric 

ds 2 = —e«dr 2 —eV (dB 2 +sm 2 ddy 2 ) + e 1 c 2 dt 2 , ( 1 ) 

where a, p, and y are functions of the radial coordinate r and the time coordinate x = ct, 
and the equation for the massless scalar field will be 

f - -1 (a' - 2p' - y') f - e'-'* 

We will assume that 

a = a 0 + 3a, p = p o + Bp, T = 7n +o T , ,i = fc, + (3) 

where a 0 » g 0 > Yo > and i|) 0 are the solution of the Einstein equations for static fields of 
a point source, and 6a, 6p, , 6y» and 6 are small perturbations that are functions of time. 
We will obtain the equation for 6ij), leaving the form of the static solution indeterminate 
(in contrast to [1]). We will select the perturbation coordinates so that 6(3 = 0. Sub- 
stituting (3) into (2) gives the equation for 6^, which contains the quantity (6a 1 — 6 > ' ) . 
We will express this quantity in terms of 6(p. We obtain the following from Einstein's 
equations for metric (1) with the energy-momentum tensor of a massless scalar field (as- 
suming that (J = Pq and is not a function of time) 
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T)<H=0. 



(2) 



and 



~P'K-7') = r + 3' a -2e«-3. (5) 

Describing the latter equality first for the static solution and then for the "perturbed" 
solution while using (4), we can obtain the quantity (6a 1 — <5y'), after which we will find 
the equation for a small perturbation in the potential of a massless scalar field in the 
form 

W ~ j f (r) 3f + f(r) W - g (r) 8 $, (6) 



where 



(7) 



Separating variables, we have 



g (r) = e°.-T., ( 8 ) 

?(r) = a;-2p;- T ;. (9) 

= 8\to(JC°)6\|)i(r), 6%, = exp(±(W), (10) 
8-K - I T (r) 8ft + (o)-^ + /) Sti = 0. (11) 

The separation constant (qj 2 ) determines the circular frequency (wc) of the perturbed 
radial "pulsations." If (11) has negative eigenvalues u 2 , then, according to (10), it is 
possible to have an unlimited increase of 6il) 0 with time, which is also true for <5i|j, i.e., 
the scalar field in this case is unstable. 

The actual form of Eq. (11) depends on the form of the static solution. The form of 
the latter is determined by the selection of the coordinate conditions for which the static 
Einstein equations are solved. 

We will select the following relation for the coordinate condition 

do— 2p 0 — To = — 41nr. (12) 

Then, the solution of the static equations in the Newtonian limit will be 

p 0 = .(*+ 1)^ + 2 In [*r f / 

1>o == - q\r. To = - r g (r, 

where r g is the gravitational radius of the source, q is its scalar charge, and 



(13) 



In this case 



/(/■)== g(r)= —, «p(r) = -±. (15) 



Introducing the new function 
we obtain the equation 



y{r) = rtyu 



Assigning homogeneous boundary conditions at the origin and at infinity, we will obtain an 
eigenvalue problem whose solution answers the question of the field's gravitational sta- 
bility. 
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Equation (16) was obtained by the methods indicated in [1], where the stability of the 
static solutions of the Einstein equations is investigated for small perturbations. However, 
it differs from the corresponding equation in [1] because a different form of the static 
solution was used (it can be obtained from the Bronnikov solution [2] through a coordinate 
transformation). One should note that the Bronnikov solution has many advantages over the 
Fisher solution because the components of the metric tensor and the potential of the scalar 
field are expressed by direct functions of an independent variable). 

The eigenvalue problem is also formulated differently. Since Eq. (16) possesses no 
singularities, the boundary conditions are first order, which requires that the function 
y(r) go to zero at the origin and at infinity. In [1], the equation was first transformed 
into the normal Liouville form. This transformation of Eq. (16) is completely unnecessary 
for our case, since it "degrades" the equation, giving rise to singularities and making it 
difficult to formulate single-valued boundary conditions in simple form. This contrast in 
the formulation of the eigenvalue problem also leads to results in our study and in [1] 
that differ significantly. 

We will estimate the boundaries in which the quantity of the smallest eigenvalue to 2 
lies using the inclusion theorem in [3]. Using homogeneous boundary conditions for Eq. (16), 
the inclusion theorem confirms that the smallest eigenvalue lies between the minimum and 
the maximum of the auxiliary function 

O(r)^, (17) 
gu 

where u is an arbitrary comparison function (positive over the entire infinite interval). 
Using the dimensionless variable 

x=r/r g (18) 
we can preserve the form of Eq. (16), but instead of the quantities to and f, we have 

Q = o)-r g , ) ( 19 ) 
h=f-r g . f 

Using the static solution (13) gives the equations 

, / K fx V r 2(« + i) 

U:-f) - cxp [- 



X 



(20) 



h(x) =2k'(k 2 - I)- — . (21) 

JC 4 \(K—l)e"* + K+l]' 

From here, it is evident that at the origin (for x = 0), both g and h are equal to zero. 
When x = °°, the function h(x) is again equal to zero, but g(x) is equal to unity. Neither 
functions possess any singularities. 

We will select the following functions for the comparison functions 

' "i = TT-,> M a = — exp( ). (22) 

1+jc 2 x \ xj 

The function $ x (x), which is obtained after substituting into (17) u 1 and expressions (20) 
and (21), determines the lower limit of the smallest eigenvalue, and the function corres- 
ponds to $ x (x) and gives the upper limit. Values for the functions $ x (x), 4> 2 (x), and h(x) 
were tabulated for different values of the parameter k using a computer. It was found that 
when k changes from 1 + 10" 21 * to 2.0, the lower limit of the smallest eigen value varies 
from -0.4806 to -0.5514, and the upper limit varies from 0.2949 to 0.9445. For large values 
of the parameter k, the lower limit remains negative, and the upper limit is positive. 
Hence, although it gives a representation of the order of the smallest eigenvalue, the in- 
clusion theorem does not allow one to conclude whether the eigenvalues for different values 
of the parameter k are negative or positive. One must use the Rayleigh principle for ob- 
taining the answer to this question. In our case, the Rayleigh ratio will be [3] 
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§u"dx~ I(k) oo 
X [«] = i ■ ../(*) = J h (*) « 2 ^- 



The functions 



1 



V> + !)(» + 2) 



« 2 ^(*), 



where L£(x) are the generalized Chebyshev— Laguerre polynomials , form an orthonormal system of func- 
tions over an infinite interval (0, «) [4] . Assuming that the estimate for the first eigenvalue 
using the Rayleigh principle coincides with the null approximation of the Galerkin method, 
where the latter gives more reliable results when using orthonormal functions, we will 
select the first function of the indicated system as the comparison function in the Rayleigh 
principle, i.e. the function 

The first integral in the numerator of the Rayleigh ratio can then be calculated analytically 
fairly easily. It is equal to 1/4. The second integral (i) for different values of the 
parameter k was determined by the Gaussian method on a computer with an accuracy of 10" 6 . 
For a change in k from 1 + 10" 6 to 100, the quantity I varies from 0.932 to 1.2- 10" 6 . 

The decrease in the quantity of the integral I with an increase in the parameter < is 
due to the character of the dependence of the function h on <. Calculations shown that when 
k is close to unity, the function h(x) has a narrow, high peak near the origin. When k is 
increased, this peak is displaced in the direction of increasing x, where it becomes broader 
and its height is less. For large k, the peak disappears, and the function h(x) is nearly 
zero over the entire infinite interval, since the integral I is also small for other test 
functions. It was found that for small k, the Rayleigh ratio is negative, while for large 
values it is positive. Therefore, negative eigenvalues correspond to small k, and positive 
eigenvalues correspond to large values. The quantity k, for which the eigenvalues changed 
sign, i.e., a transition from instability to stability, was only found for the test function 
indicated above. It was established that this transition occurs when k ^ 4, where the numer- 
ator of the Rayleigh ratio is equal to zero (I = 0.25). 

It was found, then, that there is a dimensionless number which limits the region of 
stability for a massless scalar field. The role of this number in the problem of the gravi- 
tational stability of a massless scalar field is similar to, for example, the Reynolds num- 
ber in hydrodynamics. 

The instability of a massless scalar field of a point source when the ratio of the charge 
of the source to its mass is small indicates that a scalar charge that is too small cannot 
destroy the horizon in the metric, since the field is unstable and disappears. 
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